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ABSTRACT 

We present the first VLT near-IR observations of a gravitationally lensed quasar, using adaptive optics and laser guide star. These 
observations can be considered as a test bench for future systematic observations of lensed quasars with adaptive optics, even when 
bright natural guide stars are not available in the nearby field. With only 14 minutes of observing time, we derived very accurate 
astrometry of the quasar images and of the tensing galaxy, with 0.05" spatial resolution, comparable to the Hubble Space Telescope 
(HST). In combination with deep VLT optical spectra of the quasar images, we use our adaptive optics images to constrain simple 
models for the mass distribution of the lensing galaxy. The latter is almost circular and does not need any strong external shear to 
fit the data. The time delay predicted for SDSS J0806+2006, assuming a singular isothermal ellipsoid model and the concordance 
cosmology, is Af ^ 50 days. Our optical spectra indicate a flux ratio between the quasar images of A/B=1.3 in the continuum and 
A/B=2.2 in both the Mgii and in the Cm] broad emission lines. This suggests that microlensing afi"ects the continuum emission. 
However, the constant ratio between the two emission lines indicates that the broad emission line region is not microlensed. Finally, 
we see no evidence of reddening by dust in the lensing galaxy. 
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1. Introduction 

Strong gravitational lensing of distant quasars is a useful tool 
both for cosmology and for studying luminous and dark matter 
in massive galaxies. For a fixed cosmology, the total projected 
distribution of matter in the galaxy(ies) responsible for the lens- 
ing effect can be reconstructed from (i) the configuration of the 
lensed quasar images on the plane of the sky, (ii) the observed 
magnification ratio between the quasar images, (iii) measure- 
ment of the time delay, and (iv) detailed structures in the lensed 
host galaxy of the quasar forming arcs and counter-arcs (e.g. 
Refsdal 1964 Claeskens & Surdej 2002, Kochanek et al. [2^ . 
Finally, the magnification ratios between the quasar images can 
significantly deviate from the predictions by smooth mass mod- 
els if the lensing galaxy contains small-scale substructures with 
M ~ 10'*-'' Mq (e.g., Mao & Schneider 1998, Metcalf [2005]l. 
Strong gravitational lensing of quasars is therefore a sensitive 
indirect way of detecting small satellites of lensing galaxies. 

High-resolution imaging of lensed quasars is crucial for con- 
straining mass models for the lensing galaxy, determining its lu- 
minous mass distiibution, and identifying the optical counterpart 
of any mass substructure (Schechter & Moore 1993 Koopmans 
& Treu ,2002 McKean et al. ,2007j . The Hubble Space Telescope 
(HST) has provided the community with a large sample of 
exquisite high-resolution images of gravitational lenses with the 
CASTLES survey (CfA-Arizona Space Telescope LEns Survey; 
Munoz et al. 119981) and played a crucial role in the study of 
lensed quasars (Claeskens et al. 120071 for a review). Ground- 
based adaptive optics (AO) imaging now allows us to obtain 
high spatial resolution data from the ground. Although both ap- 



proaches are complementary in depth and field of view, the AO 
imaging data collected so far of lensed quasars are still very 
scarce. With Gemini-North, Courbin et al. (I2002I I obtained the 
first AO images of a gravitational lens, for the Einstein ring 
PKS 1830-21. These observations were followed up at the VLT 
in AO in order to confirm a double lensing galaxy (Meylan et 
al. ^005 ). More recently. Auger et al. (120081 ) used AO to ob- 
serve the doubly imaged quasar SBS 1520-1-53. These two tar- 
gets are easy for AO, as they are located a few arcseconds away 
from a blight natural guide star. Systematic AO imaging of all 
other lensed quasars requires laser guide star (LGS). McKean et 
al. (12007b obtained the first AO/LGS images at the Keck of a 
lensed quasar to search for substructures in the lensing galaxy. 

In the present paper, we describe the first AO/LGS images 
of a gravitational lens at the VLT: the doubly imaged quasar 
SDSS J0806+2006, at z, = 1.540. This quasar has been found 
to be gravitationally lensed from a systematic search for strong 
lens systems in the Sloan Digital Sky Survey (Inada et al. l2006l l. 
The angular separation between the two quasar images is AO ~ 
1.5" and the early-type lensing galaxy has been found to be at 
Zi = 0.573 from deep VLT optical spectroscopy (Eigenbrod et 
al. l2007l l. Our new VLT AO/LGS images, further improved with 
image deconvolution, allow us to derive milli-arcsecond astrom- 
etry of the quasar images, as well as improved models for the 
lensing galaxy. We see these observations as a test bench for 
systematic VLT AO imaging of gravitationally lensed quasars. 

2. Observations and reduction 

2. 1 . VLT adaptive optics observations 



* Based on observations collected at the European Southern 
Observatory, Paranal, Chile, ESO Program 079.A-0588(A). 



Our AO/LGS H-band images of SDSS J0806+2006 were ob- 
tained in service mode at the ESO Paranal Observatory with 
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Fig. 1. Left: VLT image of SDSS J0806+2006 obtained with the adaptive optics + laser guide star facihty. The PSF FWHM in this 
14-minute //-band image is about the same as on a HST/NICMOS image in the same band. Right: deconvolved image, as obtained 
with the MCS algorithm (see text). The lower cut level in this deconvolved image is set to Icr^i^y in order to show the noise level. 
North is to the top, east to the left. 



the near-infrared camera CONICA (the COude Near-Infrared 
CAmera), which is mounted on the AO system NAOS (Nasmyth 
Adaptive Optics System), installed at the Nasmyth B focus of the 
VLT-UT4. This instrument uses an LGS Facility, in operation at 
Paranal Observatory since April 2007. The LGS is an artificial 
source replacing the natural guide star as a reference object for 
AO image correction. A natural guide star is still required to 
correct for the tip-tilt motions, which are not sensed by the LGS, 
but the natural guide star can be much fainter (12 < V < 17) 
than without the LGS. In addition, the natural guide star can be 
as far as 55" away from the science target. The S27 CONICA 
camera was used, with a pixel size of 27.053±0.019 mas, result- 
ing in a field of view of 28x28"(Amico et al. 2008). A set of 7 
exposures of 120 s each were obtained in the //-band on 2008 
January 7, at a mean airmass of sec(z) ~ 1.4, under seeing condi- 
tions in the range 0.65"-1.0". We used the star U0975_073 17639 
(V=14.2), located 44.6" from SDSS J0806+2006, as a reference 
tip-tilt star. These relatively poor observing conditions led to a 
Strehl ratio on-target of only 6%, with a final full width half at 
maximum spatial resolution of 0.16" on the coadded image. The 
point spread function (PSF) varied from frame to frame, with an 
FWHM varying from 0.13" to 0.19" and an ellipticity varying 
between 0.0 and 0.2. 

Standard near-IR reduction procedures were applied to sub- 
tract the dark and to flat-field the images using a normalized 
twilight flat field. The sky subtraction and the co-addition of 
the reduced dithered frames were performed using the xdimstm 
IRAl|^package. The sky frames were computed using the 3 im- 
ages taken nearest in time to the frame considered. The standard 
star S301-D, observed on the same night as SDSS J0806+2006, 
was reduced in a similar way to the science data. From aperture 
photometry, we derived a zeropoint H = 23.87 ± 0.10 mag/s. 
In this calculation we adopted a Paranal extinction coefficient of 
0.06 mag per unit airmass in the //-band. The large uncertainty 



on the zeropoint reflects the inaccuracy of absolute photome- 
try with AO systems due to seeing variations, PSF halo vari- 
ations and angular anisoplanatism (e.g. Esslinger & Edmunds 
1998). The error on the relative photometry between images 
A and B of SDSS J0806+2006 should be less than 0.05 mag 
(Esslinger & Edmunds , 1998 ) due to the small PSF variation over 
the 2"x2"field covered by the lens system. 



2.2. VLT spectroscopy 



A set of 4 low-resolution spectra of SDSS J0806+2006 is avail- 
able, obtained with the FORS 1 instrument mounted on the ESO- 
VLT at the Paranal Observatory (Chile). Using these data, taken 
on 2006 April 22, Eigenbrod et al. ( 120071) measured the redshift 
of the lensing galaxy. In the following, we analyze the spectra 
again in order to infer the broad emission line flux ratio of the 
quasar images and to estimate the dust content of the lensing 
galaxy. The details of the observations and reduction of the spec- 
tra can be found in Eigenbrod et al. (2006a, 2006b, 2007 ). 

All of the observations were carried out using multi-object 
spectroscopy (MOS). The object slit encompasses images A-B 
and the lensing galaxy. The other slits were placed on 4 stars, 
that were used as flux calibrators and as reference PSFs in order 
to deblend the quasar images. The wavelength range 4450 < A < 
8650 A was covered with a resolving power R- A/AA = 200 at 
the central wavelength. 

The final combined 2D spectrum of SDSS J0806-H2006 was 
spatially deconvolved using the spectral version of the MCS de- 



convolution algorithm (Courbin et al. 2000; Sect. 3.2 1, leading to 



' IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



individual ID quasar spectra, free of any mutual contamination. 
These spectra were eventually corrected from the spectrograph 
response curve and from diff'erential atmospheric extinction us- 
ing tabulated Paranal extinction curves. 
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3. Analysis 

3.1. Imaging: deconvolution and profile fitting 

We analyzed our AO/LGS images in two ways. First, we ap- 
plied the MCS image deconvolution algorithm (Magain et al. 
11998 ). As no star was available in the immediate vicinity of 
SDSS J0806+2006, we built the PSF on the relatively isolated 
quasar image A (Fig. 1). The resolution in the final deconvolved 
image was chosen to be 0.055" FWHM. The frame was decom- 
posed into an analytical channel containing the point sources and 
a numerical channel containing an image of the lensing galaxy. 
The result of the deconvolution technique is shown in the right 
panel of Fig. 1 . We measured the position of the lensing galaxy 
using Gaussian fitting on the numerical channel of the decon- 
volution. The astrometry and photometry of the quasar images 
were returned as analytical parameters of the point-source chan- 
nel of the deconvolution. 

Second, we also analyzed the images by modeling the quasar 
images and the galaxy light distributions using GALFIT (v2.0.3c; 
Peng et al. 12002b . For this purpose we used quasar image A as 
our PSF model and simultaneously fit the 2 quasar images and 
an analytical galaxy profile, convolved with the PSF. The pa- 
rameters of the lensing galaxy are displayed in Table [T] for de 
Vaucouleurs (GdV) and exponential-disk (Exp) light profiles. 
Both models lead to equally good fits to the data, given our 
signal-to-noise, with a slightly better formal fit for the exponen- 
tial profile. However, as the spectrum of the lensing galaxy is 
typical of an early type galaxy (Eigenbrod et al. 2007 ), we use in 
the following the morphological properties derived from the de 
Vaucouleurs model. In Table[r| we give the reduced the num- 
ber of degrees of freedom (d.o.f.), the galaxy PA (east of north), 
the axis ratio q, and the effective radius R,,. The error bars asso- 
ciated to the galaxy parameters are formal error bars of the fit 
and do not include systematic errors associated with the choice 
of surface brightness profile. These systematic eiTors are likely 
to be dominant for the effective radius, R^, and axis ratio, q, as 
well as for the magnitude of the galaxy 

The deconvolution procedure and the GALFIT minimization 
lead to photometric measurements compatible within 0.05 mag 
and astrometiic measurements compatible within 0.005". We 
summarize the astrometry and photometry of SDSS J0806+2006 
in Table[2] The photometry is in the Vega system. The astrometry 
of the point sources is given by the deconvolution process, while 
the position of the lensing galaxy is given by GALFIT. Measuring 
the position of the numerical lensing galaxy of the MCS decon- 



Table 1. Results of the fit of the lensing galaxy with GALFIT for 
2 different light profiles. 



Model 


'? 

X 


d.o.f. 


PA 


1 = 


b/a 


(") 


GdV 


L06 


16372 


78.2 ± 6.8 


0.88 


±0.02 


8.9±0.3 


Exp 


L03 


16372 


83.9 ±6.8 


0.91 


±0.01 


7.0±0.1 



Table 2. Relative astrometry and //-band photometry for the 
lensed images A and B and for the lensing galaxy G. 



ID 


Aa (") 


AS (") 


//-Magnitude 


A 


+0.000 ± 0.003 


+0.000 ± 0.003 


17.19 ± 0.14 


B 


-1.316 ±0.003 


-0.701 ± 0.003 


18.00 ± 0.14 


G 


-0.918 ±0.009 


-0.514 ±0.009 


17.56 ±0.18 
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Fig. 2. Top: our deep VLT optical spectra of components A (dot- 
ted blue) and B (dashed red) of SDSS J0806+2006. Bottom: 
spectral ratio between B and A. The straight line at B/A =0.77 
(A/B=1.3) shows the average ratio estimated for the continuum 
emission. 



volution gives consistent results. The error bars quoted in Table|2] 
are Icr error bars including the estimated systematic errors. 

We find that the separation between the quasar images is 0. 1 " 
more than the one presented in Inada et al. (2006), incompatible 
with their quoted eiTor bars. An overestimate by about 6% of 
the CONICA pixel size would allow both data measurements to 
be reconciled. However, the CONICA pixel size has been re- 
cently confirmed by new measurements from Sicardy et al. (pri- 
vate communication) and by our measurements of the relative 
astrometry of 4 stars in the field of SDSS J0806+2006. 

3.2. Spectroscopy 

The original deconvolution of the 2D spectrum of 
SDSS J0806+2006 presented by Eigenbrod et al. (l2007l l is 
aimed at measuring the redshift of the lensing galaxy. We 
optimize here the deconvolution in order to improve the signal- 
to-noise ratio in the spectra of the quasar images and minimize 
the contamination of the point-like sources by the extended 
lensing galaxy. 

The spectra of the quasar images are shown in Fig. |2] 
Although the conditions were photometric during the observa- 
tions, the spectra are not corrected from the flux loss due to the 
finite slit width. Still, the relative fluxes are coiTect since the slit 
clipping is the same for the two quasar images and since FORS 1 
uses an atmospheric refraction corrector. 

The spectral ratio between the quasar images provides a 
good diagnostic of the presence of differential microlensing and 
differential extinction (e.g. Sluse et al. I2007I I. The imprint of 
quasar broad emission lines (BEL) is clearly visible in B/A 
(Fig. |2|, which is otherwise flat over the observed wavelength 
range. This is a clear signature of differential microlensing be- 
tween the continuum emission region and the broad emission 
line region. Indeed, the Einstein radius of a microlens located in 
the lensing galaxy is about R^e - 0.0012 /z"''^ y/MjM^ pc (pro- 
jected in the source plane). Consequently, the continuum emis- 
sion region, with a size < 10"'* pc < Rf,E, is more likely mi- 
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Table 3. Results of the lens models, where angles are measured in degrees east of north. 



Model 


Re{") 


7 or e 


6y or e, At{h'' days) 


ka 


7a 


Kb 


7b 


SIS+7 


0.74 


0.017 


56.8 36.3 


0.351 


0.342 


0.841 


0.833 


SIE 


0.76 


0.042 


76.6 36.2 


0.356 


0.356 


0.850 


0.850 



crolensed than the much larger broad emission line region. The 
flux ratio in the continuum is found to heA/B{cont) - 1.3, while 
the ratio in the Mgiiand Cm] is found to he A/B(BEL) - 2.2. 
Interestingly, the BEL measurement also agrees with the flux ra- 
tio measured in the //-band, i.e., A/B{H) = 2.1. This is consis- 
tent with the need for the //-band continuum emission to arise 
from a larger region than the optical continuum emission. Like 
Inada et al. (2006), we do not detect any significant chromatic 
change of A/B in our spectra. This indicates that differential ex- 
tinction between the quasar lensed images is small. Although 
microlensing could induce chromatic variations in the quasar 
flux ratios (Wambsganss & Paczyhski 1991), the absence of such 
changes in microlensed quasars is not uncommon. 

Finally, there is an apparent change in the flux ratio A/B 
between April 2005 and April 2006. Inada et al. (2006) report 
A/B{cont) ~ 1.45 in the continuum spectral ratio and in the opti- 
cal images, while we measure A/B{cont) -1.3. This is compati- 
ble with an increase of 10% in the micro-magnification between 
April 2005 and April 2006. Systematic errors in the spectral 
extraction and in the image flux measurements, along with the 
combined effects of intrinsic variability and a time delay might 
also explain the observed difference. 



4. Lens Modeling 

Due to the mismatch between the astrometry exposed in Table |2] 
and the one in the discovery paper (Inada et al. 2006 ), we present 
new mass models of the lensing galaxy. Two standard mass mod- 
els are used: a singular isothermal sphere with an external shear 
(SIS-Hy) model and a singular isothermal ellipsoid (SIE) model. 
Both models have eight parameters: the Einstein radius Re, the 
shear y or ellipticity e and their respective position angles {9y or 
6e), the position of the lens galaxy, and the position and flux of 
the source quasar We have 8 nominal constraints (the positions 
of A, B, and G, the flux of A and B). The PA of the galaxy light 
distribution can also be used to constrain the SIE mass model. 
Indeed, statistical studies find correlations between the position 
angles but not between the axis ratios of the visible and total 
mass distributions (Keeton et al. 119971 Ferreras et al. 120081 ). In 
conclusion, the SIS-i-y model has d.o.f., while the SIE model 
has 1 d.o.f. We used the quasar and galaxy positions reported in 
Table|2]and the flux ratio B/A - 0.455 ± 0.050 measured in the 
broad emission lines. We chose this flux ratio because it is less 
affected by microlensing than the continuum flux ratio. 

We fit the data using LENSMODEL (vl.99g) (Keeton |200TT ). 
The results are summarized in Table |3] In addition to the model 
parameters, we give the local convergence k and shear y at the 
location of the lensed images as derived from the magnifica- 
tion tensor The SIS+y model fits the data perfectly, as expected 
(d.o.f.= 0). The shear is found to be very small (y =0.017). The 
reduced for the SIE model is - 0.81. This model shows 
that there is excellent agreement between the PA of the light 
and of the mass distributions. These results diff'er significantly 
from those presented in Inada et al. (2006 ). The origin of this 
discrepancy is the significantly different relative astrometry of 
SDSS J0806-H2006 obtained with our AO imaging. 



5. Conclusions 

We obtained the first VLT images of a gravitationally lensed 
quasar using AO and LGS. Using image deconvolution we de- 
rived accurate astrometry of the quasar images and of the lensing 
galaxy. This astrometry is significantly different from the one in 
Inada et al. (2006) and is not due to a systematic error in the 
pixel size of the instrument used. 

We find that a simple isothermal mass model can account 
for the new astrometry, with no need of strong external shear. 
The predicted time delay for this model. At ^ 36.3/!"' days, is 
well-suited to an accurate measurement using photometric light 
curves (see simulations by Eigenbrod et al. 2005 ). 

Our VLT optical spectra show that SDSS J0806-H2006 is 
affected neither by chromatic microlensing nor by dust extinc- 
tion. The continuum flux ratio based on these spectra agrees 
within 10% with the one derived by Inada et al. one year 
before, suggesting low-amplitude microlensing-induced varia- 
tions. SDSS J0806-(-2006 is therefore a clean lens system to use 
either for determining of Hq or for the study of the mass distribu- 
tion in the lensing galaxy, based on a time-delay measurement. 

The accuracy of the relative positions of the lensed images 
derived from our observations of SDSS J0806-I-2006 is simi- 
lar to the typical accuracy derived from the HST data of other 
lensed quasars (e.g. Impey et al. 1998, Claeskens et al. 2006|. 
This demonstrates that ground-based AO imaging with the help 
of LGS can be fully competitive with HST observations in the 
near-IR. About 90% of the lensed quasars observable with the 
VLT are either bright enough to be used as a tip-tilt reference or 
have potential tip-tilt stars with V < 17 within a 55"radius. This 
is promising for future systematic observations of gravitational 
lenses with adaptative optics and laser guide star. 
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